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We study 5=1 spin liquid states on the Kagome lattice constructed by Gutzwiller projected 
Px + ipy superconductors. Depending on the topology of the fermions, the obtained spin liquids 
can be either non-Abelian or Abelian. By calculating the modular matrices 5 and T, we conhrm 
that projected topological superconductors are non-Abelian chiral spin liquid (NACSL). The chiral 
central charge and the spin Hall conductance we obtained agrees very well with the 50(3)i held 
theory predictions. The NACSL may be stabilized by a local Hamiltonian. From a variational study 
we observe a topological phase transition from the NACSL to a Z 2 Abelian spin liquid. 


Since the discovery of fractional quantum Hall (FQH) 
effect topological order became a fundamental con¬ 
cept in condensed matter®^®. In contrast to usual long- 
range orders characterized by spontaneous symmetry 
breaking according to Landau’s paradigm, topological or¬ 
ders are characterized by the topological degeneracy of 
ground states on a manifold with genus and the fractional 
bulk excitations with a nonvanishing gap. For instance, 
the FQH liquid with ^ filling has three-fold degenerate 
ground states on a torus and its quasiparticle excitations, 
called anyons, carry | charge. The charge-1 anyons re¬ 
spect fractional statistics, i.e., the many-body wave func¬ 
tion obtains a Berry phase if one anyon adiabat- 

ically exchanges its position with another one (called 
braiding). More interestingly, the Pfaffian state proposed 
by Moore and Read® for the v = \ FQH liquid^’® sup¬ 
ports non-Abelian anyons in the vortex cores. The braid¬ 
ing of non-Abelian anyons is equivalent to operating a 
matrix on the degenerate states instead of multiplying a 
C/(l) phase factor®d0. non-Abelian topological or¬ 

ders have potential applications in quantum memory and 
quantum computation^®’^^. 

Besides FQH systems, gapped spin liquids also ex¬ 
hibit nontrivial topological orders. For example, the 
Kalmeyer-Laughlin chiral spin liquid^^ supports semionic 
anyons, and the short-range resonating valence bond 
(RVB) state on a two-dimensional (2D) non-bipartite lat¬ 
tice carries Z 2 topological order^®. In seeking of spin 
liquids in realistic microscopic models, antiferromagnets 
on the Kagome lattice have been widely studied^"^^®® 
due to its strong geometric frustrations, which is neces¬ 
sary for suppressing Neel order. From experimental side, 
promising candidates of spin liquids have been synthe¬ 
sized, such as the Herbertsmithite realizing the Kagome 
antiferromagnet^®. Recently, 5=1 antiferromagnets 
have attracted lots of interest from experimental, the¬ 
oretical and numerical sides. Several exotic spin liquid 
phases, such as U{1), Z 2 , and chiral spin liquids, have 
been proposed in 5 = 1 spin systems®®^^®. 

In this paper, we construct both Abelian and non- 


Abelian 5=1 spin liquid states on the Kagome lat¬ 
tice. These wave functions are constructed by Gutzwiller 
projection of Px + fpy-superconductors in the fermionic 
slave particle representation®^’^®. It is shown that the 
topology of the fermions determines the physical proper¬ 
ties of the spin wave function after the projection. The 
projected topological superconductors are non-Abelian 
chiral spin liquids (NACSL), which is verified by calcu¬ 
lating their modular matrices T and 5. We show that 
the NACSL has quantum spin Hall effect where the spin 
Hall conductance is quantized to On the other hand, 
the projected trivial superconductors are Abelian Z 2 spin 
liquids without spin Hall effect. We study a spin Hamil¬ 
tonian for which the spin liquid states might be stabilized 
as its ground states. By tuning the interaction parame¬ 
ters, we observe a topological phase transition from the 
non-Abelian spin liquid phase to the Abelian spin liquid 
phase. 

Construction of Gutzwiller wave functions. Following 
Ref. 27, we introduce three species of fermionic slave 
particles Ci = (cii, coi, C-ii)"^, called spinons, to rep¬ 
resent the 5 = 1 operators as S^' = :^(5f + iSf) = 

c\iCoi + cjc_ij and Sf = cl^cu - under the par¬ 
ticle number constraint -I- CQ^coi + c^nC-u = 1. In 

this way, the spin-spin interactions can be represented in 
terms of interactions between the fermions. The spin op¬ 
erators and the particle number constraint are invariant 
under local 1/(1) transformation Ci —i Cie^‘^‘, indicating 
an t/(l) gauge structure of the above fermionic represen¬ 
tation. 

We will focus on 50(3) spin-rotationally invariant 
mean-field states on the Kagome lattice (see Fig. 1(a)), 
which correspond to 5 = 1 spin liquid phases. The mean- 
field Hamiltonian for fermions takes the following form: 

+ A.^cjc, + h.c.] + 0/0,(l) 

(U) i 

where Ci = (cL^j, —c|j)^. Both Ci and Cj are SO{S) 
vectors since CjCj\vac) is a spin singlet. Noticing the 
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FIG. 1. Laughlin’s gauge invariant argument of spin Hall 
conductance. The calculation is performed on a cylinder with 
10 X 10 unit cells (300 sites), (a) the structure of Kagome 
lattice; (b) the pumping process; (c) projected topological 
superconductor (with parameters x = ^ = 1) 

has quantized spin Hall conductance ^; (d) projected trivial 
superconductor (with parameters x = ~l,^ = l:'^ = 10) has 
vanishing spin Hall conductance. 


pairing term has odd parity, the pairing symmetry on the 
Kagome lattice can be either p-wave or /-wave, depend¬ 
ing on the angular momentum of the ‘Cooper pair’. Here 
we will focus on the {px + ipy)-pairing states with vari¬ 
ational parameters Xij = Xji = ^ij = ~^ji = 

{0ij is the angle between a;-axis and the bond ij) and 
the ‘chemical potential’ A. In the weak pairing region 
—2|x| < A < 4|x|, the mean-field state is a topological 
superconductor with nontrivial Chern number^®. On the 
other hand, in the strong pairing region A < — 2|x| or 
A > 4|x|, the Chern number vanishes and the mean-held 
state is a trivial superconductor. 

To go beyond the mean-held theory, we need to con¬ 
sider the low energy 17(1) phase huctuations of the varia¬ 
tional parameters, which are Higgsed into Z 2 gauge huc¬ 
tuations owing to the fermion pairing. Alternatively, a 
RVB wave function can be obtained by a Gutzwiller pro¬ 
jection of the mean-held state, 

I RVB) = PgIBCS), (2) 

where Pq is the Gutzwiller projection operator that en¬ 
forces the onsite particle number constraints^ and |BCS) 
is the ground state of the mean-held Hamiltonian (1). In 
the following we will provide evidence that the projected 
topological superconductors are non-Abelian chiral spin 
liquids. On the contrary, the projected trivial supercon¬ 
ductors are Z 2 Abelian spin liquids. 


Response field theory and Integer Quantum Spin Hall 
effect. In the topological superconductor, each havor 
of fermions has nontrivial Chern numbers, so we ex¬ 
pect that the system has nontrivial response when it is 
probed by 5'0(3) symmetry twisting helds A^-, where 
Aij = AfjS^ + AfjS'^+AfjS^ behave like external SO{3) 
gauge helds coupling to the spins. 

Owing to the SO{3) gauge invariance, we expect that, 
after integrating out the fermions and the Z 2 gauge huc¬ 
tuations, the following SO{3)i Chern-Simons response 
theory describes the low-energy physics in the hydrody¬ 
namic limit: 

Aeff = xTr[e^'^^(A^cl,^A>,) — -A^] -I-/iMaxweii + ■•■) 

47r z o 

with k = 1. The SO{3)i Chern-Simons theory supports 
non-Abelian anyonic excitations, indicating that the spin 
wave function (2) describes a non-Abelian spin liquid. If 
the probe held only contains z component, A^ = AfS^, 
then we obtain the spin Hall current, Jf = 
where Ff = ie^''^{d^A^) is the strength of the probing 
held and the spin Hall conductance is quantized to 

We can study the response of the spin system from 
Gutzwiller projected wave functions. From Laugh- 
lin’s gauge invariant argument^^, the spin Hall con¬ 
ductance can be obtained by measuring the spin 
pump of the Gutzwiller projected state in a cylin¬ 
der when adiabatically inserting a symmetry hux 
(see Fig. 1(b)) in the mean-held Hamiltonian 

^mf(<l>.) = + 

h.c.] -I- A CjCi, where ^ A^’ • da; is the symmetry 
hux through the cylinder. Our numerical data for the 
projected ground state of H^f{^ 2 ) is shown in Fig. 1(c), 
where the spin Hall conductance is given by 
Alternatively, the spin Hall conductance can be olotained 
by calculating the many-body Chern number on a torus 
with twisted boundary conditions'^. When using the lat¬ 
ter scheme, the many-body Chern number we obtained 
is also 1, which agrees with the result from Laughlin’s 
gauge invariant argument. 

Similar to electronic quantum Hall states, the bulk of 
the NACSL is gapped and the boundary is gapless. This 
is verihed by computing the spin-spin correlation func¬ 
tions (see Supplemental Material). 

Ground-state degeneracy and modular matrices. Now 
we show that the projected topological superconductors 
are Ising-like non-Abelian spin liquids. Without lose of 
generality, in the following we will consider projected 
state with parameters x = —1,A = 1,A = 1.5. The 
Ising-like topological order contains three kinds of anyons 
(/, CT, tf) with fusion rules axa = I-\-'4>, tljX'tf = I, 
a X Ip = a. Here / is a trivial anyon which stands for 
the vacuum, ^ is a fermion and a is the Ising-like anyon 
(this topological order is named Bef. 33). 

We first figure out the ground-state degeneracy. Notic¬ 
ing that inserting Z 2 gauge fluxes through the two 
holes of the torus (which is equivalent to switching 
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the fermion boundary conditions from periodic ones to 
anti-periodic ones) does not change any local physical 
properties of the spin system, so we have 4 different 
mean-field states labeled by the Zi fluxes in the holes 
(0, 0), (0, tt), (tt, 0), (tt, tt). However, not all these states 
survive after the Gutzwiller projection. In fact, the state 
(0, 0) has odd fermion parity, so if the lattice sites are 
even, this state vanishes after projection |Pg(0,0)) = 0. 

The three remaining states survived after Gutzwiller 
projection are degenerate in energy as long as the spin 
interactions are short-ranged. Now we further show that 
these three states are orthogonal to each other. To this 
end, we calculate the overlap between these three states 
on a torus with 10 x 10 unit cells. With 3 x 10® steps of 
Monte Carlo simulation, we obtain 

|(PG(0,^)|PG(7r,^))| = 0.0012, 

|(PG(7r,0)|PG(7r,^))| = 0.001, 

|(PG(7r,0)|PG(0,7r))| = 0.0002. 

These results indicate that the three states are orthonor¬ 
mal up to errors of order 10“®. So we conclude that the 
degeneracy of ground states on a torus is three. This re¬ 
sult supports the Ising-like topological order of the sys¬ 
tem since the degree of ground-state degeneracy on a 
torus should be equal to the number of simple anyon 
types. 

To completely identify the topological order, we need 
to calculate the modular matrices S and 7^34,35 ^ 

T are the representation matrices of the corresponding 
modular transformations S and T on a torus, where S is 
a 90° rotation and T is a Dehn twist. The two transfor¬ 
mations S', T generate the modular group of a torus. We 
apply the method proposed in Ref. 36 to calculate the 
Tx and Ty matrices, namely, the representations of the 
Dehn twist along x- and j/-directions respectively. The 
basic idea is to adiabatically perform the Dehn twist with 
many sub-steps and then calculate the Berry phase in the 
whole process via wave function overlaps. To enhance the 
overlap of wave functions, each step of the Dehn twist 
is realized by shifting the boundary-crossing couplings 
by one lattice site (also see the Supplemental Material). 
Generally, the size-dependence of the Berry phase (in our 
calculation we have set = Ly = L, where and Ly 
are the number of unit cells along x- and ^-directions, 
respectively) is given by 

ip{L^)=a + ^L^ + 0{L-^). 

The intercept a is universal and only depends on the 
topological order. 

Notice that the Dehn twist can interchange or, in gen¬ 
eral, mix the degenerate states in the ground-state sub¬ 
space. For simplicity we denote \Pg{t^,0)), |PG(0,7r)) 
and \Pg{'x,'x)) as jcc), \y) and \xy) respectively. The 
Dehn twist Ty keeps the sector jx) invariant and ex¬ 
changes the sectors \y) and \xy). We denote ipx,x as the 
Berry phase of |a:) during the Dehn twist. Since the Berry 



FIG. 2. (Color online) Scaling of Berry phase for Dehn twist 
along y-direction. We have set = Ly. The meaning of the 
subscript x, y and xy in y} are explained in the main text. 

phase is not well defined for an open path, ipy^xy is de¬ 
fined as half of the Berry phase during a double Ty twist 
(namely, Ty) for the state \y). The results of the scal¬ 
ing of <p{L'^) are shown in Fig. 2. From the data we get 

Ty ={ 0 0 |, where ax^x = 0.23687r and 

V 0 0 j 

cty,xy = o.xy,y = —0.132l7r. Here we omitted the over¬ 
lap between different topological sectors. Owing to the 
equivalence between x- and y- axes, we can write out Tx 

as 7 'x — I 0 0 I, tvhere fXx.xy — ^y.xy and 

Vo 0 

^y,y ~ ~^x,x (these relations have been numerically ver¬ 
ified). The S matrix can be obtained from the relation 
S = TyT-%. 

The T and S matrices can be transformed into the 
standard form with a unitary transformation 

/I 0 0 \ / 1 V2 1 \ 

= 0 0 , 5 =— ^ 0 -V2]. 

Vo 0-1/ V 1 -^/2 1 ) 

y 1331/^ ) — g 0.02T47r2 

In this form, the diagonal entries of Ty stand for the 

self-statistics angle of anyons (the prefactor is cho¬ 

sen such that the statistics for the trivial and fermionic 
anyons are exact), and the entries of S are proportional 
to the mutual statistics angle between different anyons. 
The phase factor in S is owing to numerical errors. 
From the S matrix, we can read out that the quantum 
dimensions of the anyons are l,-\/2,l, which are consistent 
with those of /, tr and ip, respectively. 

From the Ty and S matrices, we can obtain 
the chiral central charge via the relations Ty = 
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FIG. 3. (Color online) (a) two 7r-vortices are connected with a 
string, all the mean-field couplings crossing the string reverse 
their sign; (b) correlation of the ‘cluster spin’ S 3 (i) = ’ 

which is defined as the total spin moment of the three vertices 
of a triangle. The red (blue) line shows the case where a tt- 
vortex is present (absent) at the center of each cluster. 

e~^‘^-diag{9i,9cr,9^) {9i is also called topological spin 
for the t-anyon) and with D'^ = 

The first relation tells us c_ = ^ay^xy = 1-59, 
while the second relation gives c_ = ^{ux^x — Oiy,xy) = 
1.48. This inconsistency is again owing to numerical er¬ 
rors. The averaged central charge is c_ = 1.53 ± 0.06, 
in excellent agreement with the SO{3)i conformal field 
theory prediction c_ = |. 

Non-Abelian any on carries fractional spin. Although 
the spin Hall conductance is not fractionalized, the non- 
Abelian anyons do carry fractional symmetry charge. To 
see this, we study the spin momentum of the vortex core 
of a TT-flux [see Fig. 3(a)], which traps three Majorana 
zero modes at the mean-field level. After Gutzwiller pro¬ 
jection, the three Majorana zero modes in the 7r-vortex 
correspond to a non-Abelian anyon (actually it is a com¬ 
posite anyon protected by the SO{3) symmetry, see Sup¬ 
plemental Material). 

We calculate the correlation of spin momentum S 3 for 
two triangles, where S 3 stands for the total momentum 
of the three spins on the vertices of a triangle. As shown 
in Fig. 3(b), if each triangle contains a tt- vortex, the cor¬ 
relation converges close to — in contrast, if there is no 
TT-vortex in the triangles, then the correlation converges 
to 0. This verifies that the non-Abelian anyon carries 
spin-1/2 similar to the edge states in S' = 1 Heisenberg 
chain^°’^^. 

Strictly speaking, the NACSL phase is an SO(3) sym¬ 
metry enriched topological phase^®, since the degeneracy 
according to the local spin-1/2 degree’s of freedom is pro¬ 
tected by SO(3). If we break the symmetry by a weak 
magnetic field, then the local degeneracy will be lifted 
and the composite anyon corresponding to the 7r-vortex 
reduces to the simple anyon a. 

A Gutzwiller projected trivial superconductor can be 
adiabatically connected to a nearest neighbor RVB state, 
which has 4-fold degenerate ground states on torus and 
contains Z 2 (Abelian) topological order. As shown in 
Fig. 1(d), this Z 2 spin liquid has no spin Hall effect, and 
spin moment of the anyons are not fractionalized. 

Local Hamiltonian and topological phase transition. 
Exact parent Hamiltonian has been proposed^®’"^® for 



FIG. 4. (Golor online) Tentative phase diagram. We have set 
Ji = 1 and only consider the region where is not too small 
and K is not too large. 


a non-Abelian spin-1 Pfaffian state. That Hamilto¬ 
nian contains three-body interactions and the interaction 
length is long-ranged. If the Hamiltonian is truncated to 
only contain short-range interactions, its ground state 
possibly falls in the same phase"*’^. Here we consider a 
simple local spin model on the Kagome lattice, 

H = ■ S, - A(S, • S,)"] -b 4 ^(S, X S,-) • Sfc,(3) 

(ij) A.V 

where Ji,K > 0, and i,j,k goes counterclockwise on 
each equilateral triangle (generally the three-body inter¬ 
actions on skew triangles also exist but are ignored here 
for simplicity). The term preserves SO(3) symme¬ 
try but explicitly breaks time-reversal symmetry, which 
may help to stabilize the NAGSL state (see Supplemental 
Material). 

Minimizing the energy of the projected Px + ipy su¬ 
perconductors with respect to the Hamiltonian (3), we 
obtain a tentative phase diagram with two spin liquid 
phases, as shown in Fig. 4. It is notable that the quantum 
phase transition between the two spin liquids is a topo¬ 
logical transition since the two phases have the same sym¬ 
metry and only differ by topological orders. We also use 
an Abelian chiral spin liquid state'^^’'*^ (projected Ghern 
band insulator) as a trial wave function and find that its 
variational energy is generally higher than the NAGSL. 
However, we cannot rule out the possible existence of 
symmetry breaking phases^^’^® which may appear in the 
phase diagram when K is large enough. The full phase 
diagram of the Hamiltonian (3) is an interesting open 
issue and deserves further investigations. 

The non-Abelian anyons can be used for topological 
quantum computation. A natural question is how to trap 
and control the non-Abelian anyons. Since 7r-vortices 
traps Ising-like anyons, the question becomes how the 
TT-vortices locally change the spin interactions. Noticing 
that the 7r-vortex in a triangle reverses the sign of the 
three-body interaction (see Supplemental Material), the 
non-Abelian anyons can be created locally by defect tri¬ 
angles with a —J^ interaction^®. 

The method of studying the Ising-like topological order 
in the present work can also be applied to study other 
topological orders. 
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Appendix A: Projected topological superconductor 
and its effective (response) field theory 


Since the fermions are interacting, the mean-field pa¬ 
rameters Xj ^ s-nd A are subject to fluctuations. The 
amplitude fluctuations are gapped, at low energy we only 
need to consider the phase fluctuations, which behave like 
gauge fields coupling to the fermionic spinons. Owing to 
the fermion pairing, the 17(1) gauge fluctuations are fur¬ 
ther Higgsed and only Z 2 gauge symmetry remains. So 
we only need to consider the Z 2 gauge fluctuations, 


then the mean-field ground state becomes 


|BCS) = n [1 “ ay (c^c^ -f -f C^c^)] |vac) 


i>j 


=n(i 

>s 

n(i 




.>11(1- 


^uvCyu^yv) 


r>s 

X 


u^v 


- ap94p4g)|vac) 


p>q 

which is essentially three copies of Px + ipy superconduc¬ 
tors. The projected states can also be written as 


Hraf{o'ij,6Xi) — -f AcjCj -f h.c.] 

{ij) 

+ ^(A + <5A,)(C|Q-1), (Al) 

i 


where uy are spatial components of Z 2 gauge fields and 
the temporal fluctuations 5Xi is still continuous. When 
the gauge fluctuations are integrated out, we obtain the 
low-energy effective field theory of the spin system. 

A practical way to obtain a spin wave function from 
the mean-field state is to enforce the particle-number con¬ 
straint by performing a Gutzwiller projection, 

|RVB) = PgIBCS) 

= ^detA(l,-l)pf 5(0,0)|a), (A2) 


where |a) is an Ising configuration, Pq is the Gutzwiller 
projection operator that enforces the particle number 
constraint, |BGS) is the mean-field ground state 


|BCS) = n + ab(4ic4j - cliclj + cl^Aj) |vac) 


i>j 


here Oy = —ciji is the relative wave function of the 
spinons in a Gooper pair. The matrices A(l,—1) and 
5(0,0) are given as 





^min2 

^m2ni 

^m2n2 


V 


5(0,0) 




where mi,ni,pi are the positions of ci, c_i, cq fermions in 
the configuration ja), respectively. The particle number 
of Cl is equal to that of c_i since the total spin moment 
Sx of the ground state is 0. 

If we introduce three new bases 


4 = - c\), cl = 4^(c 4 + cl), cl = cl 


42 ' 


42 ' 


0’ 


|RVB) = PgIBCS) 

= ^4,a;)pf D(p,p)pf B(z,z)la'} 

Ct' 

where a' is the spin configuration created by cl, cl, cl- 
The matrices B,C,D are defined as B{z,z) = 5(0,0) 
and 

( 0 

— amami 0 

( 0 ^711712 

— 0.712711 0 

here mi,ni,Pi are the positions of the Cx,Cy,Cz fermions 
in the configuration jo;'), respectively. The particle num¬ 
bers of Cx,Cy, Cz are arbitrary, the only constraint is that 
their sum is equal to the number of sites. 

If the superconducting bands have nontrivial Chern 
numbers, then ay oc |ri — at large distance and 

the projected spin wave function resembles the pfaSian 
states given in Ref. 22 and 28. 

In the topological superconductor, since each flavor of 
fermions have nontrivial Ghern numbers 1/2 (the half- 
quantized Ghern number is owing to particle-hole redun¬ 
dancy), we expect the system has nontrivial response 
when it is probed by symmetry twisting fields. Based 
on the mean-field theory, the spinons will couple to not 
only the gauge fluctuations but also the symmetry twist¬ 
ing fields (behaving like gauge fields) 

5^f(ay, Ay) = 5] ay[xyC|e*^-^Q -I- AyC^e*^-^G 
{i.j} 

+h.c.]+^(A, + (5A)(C|Q-l), (A3) 

i 

where Ay = Ay 5^-f A75^-f Ay 5^ is the external 50(3) 
probing field corresponding to spin rotation symmetry. 
Owing to 50(3) gauge invariance, we expect the follow¬ 
ing 50(3)i Ghern-Simons response theory (in the imag¬ 
inary time formalism and in the continuum limit) after 
integrating out the internal gauge fluctuations and the 
fermions, 

Tleff = i ~.— xTr[e^'^^(A^9[/A>,) — —A^] -l-5Maxweii + 

47r z o 
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where fc = 1 is twice of the fermion Chern-number since 
both the Cl and c_i fermions carry symmetry charge and 
contribute to the response. The S'0(3)i Chern-Simons 
theory supports Ising-like non-Abelian anyonic excita¬ 
tions, indicating that the spin ground state is a non- 
Abelian spin liquid. The 5'0(3)i Chern-Simons action 
is not gauge invariant under SO{3) gauge transforma¬ 
tion if space-time has a boundary. The gauge anomaly 
can be canceled by a 5'0(3)i chiral Wess-Zumino-Witten 
term, which can be represented as three species of chiral 
Majorana fermions at the boundary. In other words, the 
boundary is described by 5'0(3)i conformal field theory 
and carries chiral central charge c_ = |. 

If the probe field only contains z-component, = 
A^S^, then we obtain the spin Hall current, 

jz _ ^ pz 

^ ~ 5A- ~ 2tt 

H- 

where = ie^'^^{duAf^) is the strength of the probing 
field and the spin Hall conductance is 1 in unit of No¬ 
tice that the spin Hall conductance for bosonic 17(1) sym¬ 
metry protected topological phases and 5” = I Abelian 
chiral spin liquids are even integers time ^^^77-49^ 

Appendix B: Correlation function and zero modes in 
the non-Abelian chiral spin liquid 

From the discussion above, the response field theory 
predicts that the bulk of the NACSL is gapped and the 
boundary is gapless. We checked this numerically and 
the results are shown in Fig. 5. The spin-spin correlation 
in the bulk fits well to an exponential function, which 
confirms that the bulk is gapped; while the spin-spin 
correlation on the boundary fits better with a power-law 
function, in agreement with the fact that the boundary 
is gapless. 

Now we focus on the 7r-vortex. According to Ref. 29, 
the TT-vortex core traps Majorana zero modes in a topo¬ 
logical superconductor. In our case since the three copies 
of fermions Cx,Cy^Cz have the same mean-field Hamilto¬ 
nian owing to the S'0(3) symmetry, a 7r-vortex will trap 
three Majorana zero modes 'yx-,ly,lz (where {7m,7n} = 
2<5mn) in the weak pairing phase. In the following we 
show that these three Majorana zero modes vary as spin- 
1/2 under spin rotation. 

The three Majorana operators form an iS'0(3) vec¬ 
tor. For example, under the rotation the operators 

varies as 

7a;\ / cos 6 sin6> 0\ /jx\ 

ly I = I — sin 0 COS0 0 j ( 71 /) ,(B 1 ) 
Izj V 0 0 1] \jzJ 

where Sz is the operator acting on the degenerate Hilbert 
space of degenerate of the Majorana modes. It can be 
checked that the following expression satisfies (Bl): 



Similarly we have Sx = —^^yjz and Sy = —Ifizlx- Then 
we have 

F F si = ^ = -X {- + 1), 

from which we can conclude that the spin quantum num¬ 
ber of the Hilbert space for three zero modes is 1/2. 

Since the Gutzwiller projection does not change the 
quantum number, after projection the spin-1/2 degrees 
of freedom remain at the 7r-vortex core. To verify this 
result, we numerically calculated the cluster spin corre¬ 
lation function for two 7r-vortices, as shown in Fig. 3 
in the main text. The two vortex-cores can form a spin 
singlet or a spin triplet. If the distance between them is 
by far larger than the correlation length, then the singlet 
and the triplet are degenerate in energy. In Fig. 3 the 
calculation is performed in the singlet state. 

Notice that each 7r-vortex not only carry a local spin- 
1 /2 degrees of freedom, but also carry a nonlocal Hilbert 
space similar to the Ising anyon. Generally, braiding of 
the two TT-vortices results in a non-local operation (which 
generates entanglement) together with a local spin rota¬ 
tion for the spin-1/2 degrees of freedom. Actually, since 
the TT-vortex is shared by the three copies of Px + ipy 
topological superconductors it is essentially a composite 
anyon : crX(TXcr=(l-|-T/’)xcr = cr-|-cr. If the SO{A) 
symmetry is explicitly broken by a magnetic field, the 
local degeneracy will be lifted and the composite anyon 
will reduce to the simple anyon a. 


Appendix C: Calculation of the modular matrices 

The T matrix can be obtained from the universal wave 
function overlap 

where T is a Dehn twist operation of the torus. Generally, 
the overlap after a full Dehn twist is exponentially small 
with the increasing lattice sites. To solve this problem, 
we apply the trick introduced in Ref. 36 to divide the 
Dehn twist into many substeps such that in each substep 
the wave function changes adiabatically and the Berry 
phase can be obtained with a relatively high accuracy. 

The Dehn twist is first separated into Ly steps, in each 
step, the couplings crossing the x-boundary is shifted 
along y-direction by one lattice site(see Fig. 6) such that 
the original Hamiltonian i7[(La;,T), (l,i -|- t)] transforms 
into H[{Lx,i)A^A + t + !)]■ Noticing that during the 
shifting, some bonds cross the y-boundary and the sign 
of the coupling is affected by the y-boundary condition. 
After a full Dehn twist Ty, the boundary conditions 
(0,7r),(7r, tt) are shifted into (tt, tt), (0, tt) respectively, and 
the boundary condition (tt, 0) remains unchanged. 

Secondly, to make the twist more smooth, each step 
is further divided into several sub-steps by a parameter 





X 

(a)Correlation on the boundary. 



(b)Correlation in the bulk. 

FIG. 5. The correlations are calculated on a cylinder with 
Lx X Ly unit cells, where Lx = 20, Ly = 10 and a;-direction 
has periodic boundary condition while j/-direction is open, (a) 
The spin-spin correlation function on the (upper) boundary of 
the chiral spin liquid decays in power law; (b) The spin-spin 
correlation function in the bulk decays exponentially. 



FIG. 6. (Color online) One step of the Dehn twist operation 
Ty. We have deformed the Kagome lattice with Lx = Ly into 
a square. The blue dashed lines stand for the a;-boundary 
couplings before the twist, and the red solid line means the 
a;-boundary couplings after one step of twist. Some of the red 
lines cross the j/-boundary and the signs of the corresponding 
couplings are modulated by the y-boundary condition. As a 
result, after a full Dehn twist, the x-boundary condition will 
be modulated by the j/-boundary condition. 

However, for the sectors (tTjTt) and (0,7r), since a full 
Dehn twist exchanges them, the path is not closed. If we 
act the full Dehn twist twice, the path will be closed and 
the Berry phase will be well define. The Berry phase for 
half of the loop (a single Dehn twist) is defined as half of 
that of a double twist. 

Finally, the Dehn wist Tj, can be transformed into T_y 
(the inverse of Ty) by a global 90° clockwise rotation. The 
rotation does not affect the Berry phase, so we obtain the 
following relations: 


A G [0,1]: 

H{t,X) = H[[Lx, i),(l,i-l-t)] + XH[[Lx, i),(l,i-|-t-l-l)]. 

(Cl) 

In our numerical calculation, we take discrete values A = 
0,0.25,0.5,0.75,1. 

The Berry phase is obtained from the overlap of the 
projected ground states of above Hamiltonians 


(/)(t. A) = Im \n ({t, XlP^Pclt, X + SX 


where \t, A) is the mean-field ground state of the Hamil¬ 
tonian H (t. A) and the total Berry phase is given by 

i,A 


For the sector (tt, 0), after a full Dehn twist Ty defined 
above, the Hamiltonian goes back to itself. So the evo¬ 
lution path is closed and the Berry phase is well defined. 


Ty^y — Tx^xy — Ty,xy- 


Appendix D: Decoupling the spin-spin interactions 

In this section we explain why we can use the mean- 
field parameters x, A as variational parameters. It is 
known that the two-body spin interactions can be written 
in forms of fermions as^°’®°: 

Si • Sj = —{xljXij + 

(s, . s,)2 = (Di) 

where Xij = the mean-field 

decoupling of the Ji and K terms naturally give rise to 
the mean-field parameters x iXij) A ^ (Ay). 

Now we focus on the three-body interactions (Si x S_,) • 
Sfc. In the Cx,Cy, Cz bases, the hopping and pairing oper- 
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ators are given as 


Xij — 


^ij — ^ai^aj ; 


and the spin operators can be written as 


The three-body term can be written as 

(S, X S,) • Sfc = 5 ^ 

— ■^Xa J t ^ .„t 


(D2) 


Using the relations 




S° = -ie' 


‘/ 37 U 


CpCj, 


where a,l3,j = x,y,z and repeated indices are summed 
(the same below). 


and 

^jSxi^^p'y^TjtJ ^pa^7}-y^ 
(^P7^77A ^pA*^?77) 

Eq. (D2) can be simplihed as 


(S, X Sp • Sfc = ' 


,XijXjf.Aki 


+ AijXjkA^ki ^ijAjkXki XijXjkXk^ h.C. 


(D3) 


So using the two parameters x and A we can decouple the 
three-body interaction J^i^i ^ Sj) • in a spin liquid 
phase. Since the three-body term breaks time reversal 
symmetry, it favors a time reversal symmetry breaking 
state. If the ground state does not break spin rotation 
symmetry, then the projected Px + ipy superconductor 
will be its candidate ground state. 

Furthermore, from the expression (D3), if there is a 
TT-vortex in the triangle (ijk), then the values of (xij) 


and (Ay) on one of the three bonds (namely, ij, jk and 
ki) reverse its sign. As a result, the value of the three- 
body interaction term on the triangle {J-^(Si x Sj) ■ S^} 
also reverses its sign and the state is in a higher en¬ 
ergy comparing with the ground state. If we reverse the 
sign of on the triangles which contain 7r-vortices, then 
{—J^(Si X Sj) ■ Sk) will still have a low energy and the 
state with anyons localized near the vortices becomes the 
ground state of the new Hamiltonian. 






